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ABSTRACT 


The  ASR-9  Wind  Shear  Processor  (WSP)  is  intended  as  an  economical  alternative 
for  those  airports  that  have  not  been  slated  to  receive  a  Terminal  Doppler  Weather  Radar 
(TDWR),  but  have,  or  will  be  receiving,  an  ASR-9  radar.  Lincoln  Laboratory  has 
developed  a  prototype  ASR-9  WSP  system  which  has  been  demonstrated  during  the 
summer  months  of  the  past  3  years  in  Orlando,  Florida.  During  the  operational  test 
period,  microburst  and  gust-front  warnings,  as  well  as  storm  motion  indications,  were 
provided  to  the  Air  Traffic  Control  Tower  in  real-time. 

The  ASR-9  Microburst  Detection  Algorithm  (AMDA)  is  based  on  the  earlier 
TDWR  Microburst  Detection  Algorithm,  but  has  been  substantially  modified  to  better 
match  the  particular  strengths  and  weaknesses  of  the  ASR-9  rapid-scanning  fan-beam 
radar.  The  most  significant  additions  include  a  capability  to  detect  overhead 
microbursts,  a  reflectivity  processing  step  used  to  help  detect  velocity  signatures  that 
have  been  biased  by  overhanging  precipitation,  and  a  modification  to  some  of  the  shear 
segment  grouping  and  thresholding  parameters  to  make  them  a  function  of  range.  This 
is  necessary  to  better  accommodate  the  ASR-9.  In  addition,  AMDA  has  been  designed 
to  be  as  efficient  as  possible  to  allow  it  to  run  at  the  radar’s  4.8  seconds/scan  antenna 
rotation  rate  on  a  single  board  computer.  A  detailed  description  of  AMDA,  as  well  as  the 
performance  evaluation  strategy  and  results,  are  presented  in  this  report. 
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1.  INTRODUCTION 


The  ASR-9  Wind  Shear  Processor  (WSP)  is  being  developed  as  an  economical  alternative  to  the 
Tei  niinal  Doppler  Weather  Radar  (TDWR)  tc>r  those  airports  not  slated  to  receive  a  TDWR  system. 
Lincoln  Laboratory  has  been  sponsored  by  the  FAA  to  develop  a  prototype  version  of  the  WSP,  in¬ 
cluding  both  the  hardware  and  the  wind  sheat  detection  algorithms.  Lincoln  Laboratory  s  ASR-9 
testbed  has  been  used  to  collect  data  in  Huntsville,  AL  {19H7-8K),  Kansas  City,  KS  (1989),  and  Or¬ 
lando,  FL  ( 1 990-92),  and  these  data  have  been  used  todevclop  and  validate  the  wind  shear  detection 
algorithms.  In  addition,  during  the  summer  months  of  1990, 1991,  and  1992,  theprototype  WSP  was 
used  to  conduct  operational  demonstrations,  during  which  microburst  and  gust-front  wtumings,  as 
well  as  storm  motion  indications,  were  provided  to  the  Orlando  Air  Traffic  Control  Tower  in  real 
time. 

The  ASR-9  Microburst  Detection  Algorithm  (AMDA)  is  based  on  the  earlier  TDWR  Micro¬ 
burst  Detection  Algorithm,  but  has  been  modified  to  be  better  suited  to  the  strengths  and  weaknesses 
of  the  ASR-9  radar.  In  particular,  the  rapid  scan  rate  of  the  ASR-9  can  provide  a  faster  response  tmie 
to  microburst  events,  while  the  wide  elevation  beam  of  the  ASR-9  as  compared  with  the  TDWR 
makes  it  more  difficult  to  detect  shallow  microburst  outflows,  and  can  sometimes  result  in  enoneous 
velocity  signatures  due  to  winds  aloL.  The  typical  on-airpon  siting  of  the  ASR-9  has  also  provided 
some  new  challenges,  such  as  the  detection  of  overhead  microbursts  and  the  detection  of  events  at 
short  range  that  span  a  larger  azimuthal  extent  than  those  detected  by  an  off-airport  TDWR.  AMDA 
also  contains  new  logic  to  address  several  known  deficiencies  in  the  TDWR  algorithm.  For  example, 
the  TDWR  algorithm  can  sometimes  detect  excessively  long  shear  segments  because  the  “slope”  of 
the  shear  segment  is  only  being  validated  in  the  vicinity  of  the  segment  endpoints.  AMDA  overcom¬ 
es  this  problem  by  utilizing  a  sliding  window  to  validtue  the  slope  along  the  entire  segment  length. 

This  report  describes  the  ASR-9  microburst  detection  requirements  and  the  processing  steps  that 
make  up  the  ASR-9  Microburst  Detection  Algorithm.  Following  the  algorithm  description,  a  sepa¬ 
rate  chapter  describes  the  performance  evaluation  process  and  results.  A  final  section  discusses  the 
current  status  of  the  algorithm  as  well  as  its  future  development. 


2.  ASR-9  MICROBURST  SIGNATURE  CHARACTERISTICS 


A  microburst  is  a  sudden  downward  burst  of  wind  which,  upon  impact  with  the  ground,  results  in 
a  divergent  wind  pattern.  This  situation  is  depicted  in  Figure  1 .  As  can  be  seen  from  the  figure,  de¬ 
tecting  the  low-levc!  wind  shear  with  a  t  ASR-^^'  radar  is  made  more  difficult  than  in  the  ease  of  a 
TDWR  by  the  cosecant-squared  elevation  beam  pattern,  which  results  in  low  level  wind  estimates 
that  can  be  biased  by  the  winds  at  higher  altitudes.  To  con'ect  for  this  situation,  a  ”dual-bcam”  ap¬ 
proach  is  used,  whereby  the  vekx:ity  information  from  the  high  beam  is  used  to  estimate  the  value  of 
the  higher-altitude  wind  bias,  and  concct  for  it  (Weber,  1989). 


The  low  level  winds  originating  from  a  microburst  produce  a  run  of  radial  velocities  that  are  gen¬ 
erally  increasing  with  radar  range.  AMDA  detects  the  microburst  using  a  two-stage  approach  in¬ 
volving  the  detection  of  the  microburst  signature  across  range,  followed  by  the  grouping  of  the  de¬ 
tected  areas  across  azimuth.  The  next  two  sections  describe  the  signatures  that  are  typically  encoun¬ 
tered  in  each  of  the  two  detection  stages,  including  those  that  can  result  in  false  alarms. 

2,1  ONE-DIMENSIONAL  DIVERGENCE  SIGNATURES 

To  be  considered  hazardous,  the  velocity  increase,  or  “shear,”  produced  by  a  microburst  must 
exceed  10  m/s  over  a  distance  of  not  more  than  4  km.  An  example  of  an  ideal  signature  is  shown  in 
Figure  2a.  The  velocity  increase  is  smooth  and  almost  linear,  and  the  endpoints  of  the  shear  “seg¬ 
ment”  are  local  minima  or  maxima. 
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Range  Range 


Figure  2.  Microburst  Shear  Signature  Variations. 


Many  signatures  very  similar  to  this  exist  in  the  actual  radar  data,  but  there  are  several  variations 
that  occur  quite  frequently  in  wind  shear  regions,  and  these  also  need  to  be  detected.  The  three  most 
common  ones  are  shown  in  Figures  2b,  2c,  and  2d.  Note  that  the  sample  segments  shown  in  the  figure 
are  representative  of  velocity  radials  that  have  already  been  spatially  smoothed.  Figure  2b  exhibits  a 
region  of  extremely  shallow  or  even  slightly  decreasing  slope.  This  signature  occurs  most  frequently 
in  relatively  we'ik  shear  segments  ( <  12  m/s  ),  which  are  present  in  both  weak  microburst  events  and 
near  the  edges  of  stronger  events.  Figure  2c  is  an  example  of  shallow  slope  near  the  segment  end- 
point(s).  Because  there  is  no  clear  cut  local  maximum  or  minimum  value  in  these  cases,  it  is  some¬ 
times  difficult  to  judge  where  the  starting  or  ending  point  of  the  segment  should  be  placed.  The  re¬ 
gion  of  shallow  slope  needs  to  be  trimmed  off  if  detected  to  avoid  an  excessively  long  shear  segment. 

The  signature  shown  in  Figure  2d  contains  a  sudden  jump  in  the  velocity  field.  Even  after  spatial 
smoothing,  some  small  discontinuities  due  to  noise  or  clutter  remain  and  must  be  handled.  The  typi¬ 
cal  jump  is  2-3  range  gates  wide  with  a  velocity  difference  of  about  +/-  3  m/s  at  its  peak.  Regions 
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with  shaiper,  single-point  spikes  arc  also  present,  but  in  regions  of  moderate  to  high  radar  reflectiv¬ 
ity  (25-65  dBz),  these  are  not  as  prevalent  as  the  smaller  jumps. 

In  regions  of  low  reflectivity  (0-15  dBz),  the  ASR-9  vckx’ity  estimate  tends  to  be  somewhat 
noisy,  and  some  of  these  regions  can  exhibit  an  increasing  velocity  trend  that  mimics  the  overall 
shear  and  length  characteristics  of  a  tme  shear  .segment,  potentially  causing  false-alami  problems. 
Twoexamples  of  this  kind  of  behavior  are  shown  in  Figures  2e  and  2f.  The  first  case  (2c)  represents  a 
false  segment  created  by  noise  spikes.  The  majority  of  the  overall  shear  in  this  case  can  be  attribuied 
to  the  shear  at  the  segment  endpoi  nts.  The  second  case  (2f)  exhibits  a  generally  increasing  trend  over 
a  suitable  range,  but  the  point-to-point  velocity  differences  are  fai‘  from  smooth. 

2.2  TWO-DIMENSIONAL  DIVERGENCE  SIGNATURES 

A  two-dimensional  representation  of  an  idealized  microbursl,  along  with  the  shear  segments  de¬ 
tected,  is  shown  in  Figure  3a.  It  is  cylindrically  symmetric  and  the  maximum  velocity  differential  is 
aligned  with  the  microburst  center.  Microbursts  do  sometimes  fit  this  model;  however,  more  often 
than  not  they  are  more  inegular  in  shape  and  strength.  A  study  performed  using  data  from  Denver, 
CO,  Kansas  City,  KS,  and  Orlando,  FL,  determined  that  the  strength  asymmetry  ratio,  the  ratio  of  the 
maximum  velocity  differential  to  the  minimum  velocity  differential  when  viewing  a  microburst 
from  all  angles,  was  on  average  1.7:1,  and  can  range  as  high  as  4.0: 1  [Hallowell,  1 993].  In  addition, 
the  region  of  maximum  divergence  is  not  always  aligned  with  the  center  of  the  microburst,  but 
instead  can  lie  off  to  one  side.  Given  such  variation,  a  microburst  detection  algorithm  does  not  have 
the  luxury  of  requiring  microburst  events  to  exhibit  a  certain  two-dimensional  velocity  profile,  but 
must  instead  be  sensitive  to  a  wide  variety  of  profiles. 

The  overall  size  of  a  microburst  also  vai  ies  significantly  over  its  lifetime,  tending  to  start  out 
small  and  expand  with  time.  The  typical  area  for  a  “young”  microburst  as  seen  by  the  ASR-9  is  be¬ 
tween  1.0-1. 5  km^ ,  increasing  to  over  2.0  km^  as  the  microburst  matures,  or  in  the  case  of  a  line 
microburst,  sometimes  exceeding  10  km^.  The  interval  between  the  initial  onset  of  a  microburst  and 
the  point  at  which  it  reaches  hazardous  intensity  can  vaiy  from  as  brief  as  30  seconds  to  upwards  of 
several  minutes.  Strong  microbursts  tend  to  reach  hazardous  intensities  quite  rapidly  as  a  direct  re¬ 
sult  of  the  higher  downdraft  velocities  within  the  storm  cell.  In  order  to  detect  these  events  in  a  timely 
manner,  it  is  necessary  to  keep  the  algorithm  latency  to  a  minimum. 

Multiple  microbursts  can  exist  in  the  same  cloud,  and  as  shown  in  Figure  3c,  they  can  be  closely 
spaced.  In  such  a  case,  it  is  important  to  resolve  them  as  separate  events  in  order  to  make  an  accurate 
estimate  of  the  microburst’s  strength.  This  is  somewhat  difficult  in  practice,  because  sometimes 
what  begins  as  a  single  microburst  event  can  evolve  into  multiple  events,  and  the  point  at  which  to 
separate  the  two  events  is  somewhat  ill-defined.  It  is  also  important  to  avoid  oscillation  between  a 
single-event  detection  and  a  multiple-event  detection,  as  this  can  result  in  inconsistent  microburst 
warnings. 

Ground  clutter  and  noise  in  the  velocity  field  also  impact  the  segment  association  process.  Figure 
3d  illustrates  a  case  of  a  true  microburst  with  multiple  missing  segments  due  to  ground  clutter  break- 
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a)  Ideal  Signature  b)  Asymmetrical  Signature 


c)  Closely  Spaced  Multiple  Microbursts  d)  Missing  Segments 


e)  Noisy  Segment  Field 


Figure  3.  Microbursi  Segment  Field  Variations. 


through.  Figure  3e  depicts  an  area  containing  a  sparse  number  of  segments  resulting  from  noise  in  the 
velocity  field.  To  be  successful,  the  detection  algorithm  must  be  robust  enough  to  allow  a  substantial 
number  of  missing  radials,  while  at  the  same  time  filtering  out  cases  where  the  overall  density  of 
detected  segments  is  low. 
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2.3  MICROBURST  DETECTION  REQUIREMENTS 


An  effective  ASR-9  wind  shear  detection  algorithm  must  tx;  able  to  reliably  identify  the  variety 
of  microburst  radial  vekxity  patterns  described  above,  while  at  the  same  time  minimizing  the  num¬ 
ber  of  false  signature  detections.  In  addition,  the  4.8  seconds/scan  update  rate  of  the  ASR-9  radar 
requires  an  efficient  algorithm  in  order  to  process  the  data  at  real-time  rates  with  the  cuirent  system 
configuration.  The  ASR-9  Microburst  Detection  Algorithm  described  in  the  following  sections  per¬ 
forms  well  with  respect  to  all  three  of  these  criteria. 
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3.  WIND  SHEAR  PROCESSOR  SYSTEM  OVERVIEW 


The  prototype  WSP  data  processing  hardware,  as  well  as  the  signal  processing  steps  pcitbmied 
prior  to  the  microburst  detection  algorithm,  have  both  had  a  significant  impact  on  the  design  and 
implementation  of  AMDA.  This  section  provides  a  brief  deseription  of  the  processor  in  order  to  clm- 
ify  some  of  the  later  discussion. 

A  block  diagram  of  the  prototype  WSP  is  shown  in  Figure  4.  The  raw  time-series  data  from  the 
radar  A/D  converters  enter  the  pr(K'cs.sor  via  a  custom  Lincoln  Lab  Interface  caid,  arc  split  up  into  6 
range  “rings,”  and  are  sent  to  the  six  Mercury  Computer  airay  processor  (MC860  AP)  boards.  The 
ba.seline  signal  processing  sequence  used  during  the  Orlando  operational  tests  consisted  of: 

1 .  Adaptive  clutter  filtering  (Weber,  1987], 

2.  Low  and  high  beam  reflectivity  calculation, 

3.  Dual-beam  velocity  calculation  [Weber,  1989], 

4.  Regressive  temporal  smoothing  filter  on  dual-beam  velocity 
product,  and 

5.  9-point  nearest  neighbor  median  spatial  filter  ( low  and  high  beam  reflectiv¬ 
ity  and  dual-beam  velocity  products). 


Workstation  Workstation 


Figure  4.  ASR-9  Wind  Shear  Processor. 


Following  these  operations,  the  dual-beam  velocity  and  reflectivity  data  are  sent  to  two  Mirar 
68030  Single  Board  Computers  (marked  “DISPLAY  CPU”  and  “AMDA  CPU”  above).  The  display 
CPU  drives  two  rapid-update  analysis  displays  which  provide  site  personnel  with  a  flexible  display 
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capability  at  the  radai'  site.  The  AMDA  CPU  is  used  for  microburst  detection  prcxcssing,  the  output 
of  which  is  sent  via  Ethernet  to  the  Geographical  Situation  Displays  in  the  Air  Traffic  Control  Tower. 
Reflectivity  and  velocity  information  aie  also  output  using  the  Ethernet  connection.  Note  that  the 
other  WSP  algorithms  such  as  gust-front  and  storm  motion  are  not  cuirently  performed  in  the  V  ME 
system  shown  above,  but  are  mn  on  separate  SPARC  workstations  that  are  connected  via  Ethernet. 

The  major  points  of  interest  with  regard  to  the  design  of  AMDA  are: 

•  The  input  reflectivity  and  velocity  cUita  have  been  spatially  and  temporally 
smoothed  within  the  airay  processor. 

•  Data  are  available  to  the  algorithm  for  every  five-second  scan,  and  it  is  desir¬ 
able  to  use  all  the  data  forevci-y  scan  to  keep  the  microburst  detection  latency 
to  a  minimum. 

•  The  algorithm  is  nm  on  a  single  68030  CPU  board  rated  at  5MIPS. 

The  fact  that  the  input  data  has  been  spatially  filtered  results  in  radials  of  velocity  data  which  arc 
generally  free  from  noise  and  require  no  additional  smoothing  in  the  AMDA  processor.  This  keeps 
the  computational  requirements  of  AMDA  to  a  minimum,  allowing  data  to  be  processed  at  real-time 
rates  using  a  single  68030  CPU  board.  Th  is  is  convenient  from  both  the  softwai  e  and  hardware  stand¬ 
points  for  the  prototype  WSP. 
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4.  ASR-9  MICROBURST  DETECTION  ALGORITHM 


4.1  OVERVIEW 

AMDA  is  an  adaptation  of  the  TDWR  Microburst  Divergence  Outflow  Algorithm  (MDOA) 
(Menitt,  1991]  that  is  designed  to  run  efficiently  within  the  prototype  WSP’s  processor.  A  block 
diagram  of  AMDA  is  shown  in  Figure  5.  The  algorithm  reads  in  polar  reflectivity  and  velocity  data, 
scans  the  radar  radials  with  respect  to  range  for  divergent  outflow  patterns,  and  groups  the  detected 
regions,  teniied  "shear  segments,”  to  fomi  microburst  events.  The  first  step  performed  by  AMDA  is 
to  buffer  the  data  into  full  scans  to  allow  detection  of  overhead  events  in  a  convenient  manner.  This  is 
followed  by  the  raw  segment  detection  process  that  scans  the  individual  radials  for  microburst  signa¬ 
tures.  In  order  co  become  valid  shear  segments,  any  detected  signatures  arc  subjected  to  a  variety  of 
segment  length  and  divergent  outflow  strength  tests.  The  valid  raw  segments  are  then  temporally 
smoothed  to  provide  a  more  consistent  output  for  the  subsequent  AMDA  operations. 


POLAR 

REFLECTIVITY  & 
VELOCITY  DATA 


Figure  5.  Microburst  Detection  Algorithm  Block  Diagram. 


The  next  step,  segment  association,  groups  neighboring  segments  into  raw  “clusters,”  requiring 
that  each  raw  cluster  contain  a  certain  minimum  number  (AMDA  parameter)  of  segments.  The  raw 
clusters  are  then  passed  to  the  reflectivity  thresholding  module,  where  segments  located  primarily  in 
regions  of  low  reflectivity  (and  therefore  containing  suspect  velocity  values),  are  edited  out  of  the 
cluster.  The  cluster  is  again  required  to  possess  the  minimum  number  of  segments,  and  clusters  that 
meet  the  requirement  are  passed  to  the  segment  density  thresholding  step.  This  step  removes  regions 
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of  the  cluster  where  segment  density  is  below  a  preset  minimum,  indicative  of  regions  that  have  been 
prcxluced  by  a  few  spurious  segment  detections.  The  nKxlificd  cluster  is  again  subjected  to  the  mini¬ 
mum  segment  requirement  before  being  passed  to  the  area  thresholding  operation.  The  area  thresh¬ 
olding  operation  rejects  any  clu.sters  with  insufficient  area  to  be  considered  a  true  microburst  event, 
and  the  output  clusters  are  then  temporally  smoothed  to  increase  the  scan-to-scan  coherence  of  the 
detected  microburst  events.  An  example  of  a  "classic”  microburst,  along  with  the  velocity,  reflectiv¬ 
ity,  and  detected  shear  segment  infomiation.  is  shown  in  Figure  6.  Note  that  the  shape  shown  in  Fig¬ 
ure  6  is  for  the  rapid-update  analysis  displays  only.  Shear  segment  data  are  pas.scd  through  an  addi¬ 
tional  shape  algorithm  [Wilson,  1992]  to  produce  the  final  GSD  output. 

There  aie  a  number  of  variable  site  parameters  (e.g.  length  thresholds,  time  constants,  etc.)  asso¬ 
ciated  with  AMDA.  The  majority  of  the  parameters  in  the  following  discussion  will  be  represented 
by  symbolic  names,  followed  by  the  cunent  actual  value  where  appropriate  (e.g. 
SEG_MIN_LENGTH  (0.96km) ).  Each  of  the  parameters  has  been  optimized  using  data  from  Lin¬ 
coln  Laboratory’s  ASR-9  testbed. 

4.2  SCAN  BUFFER 

The  input  data  stream  is  in  the  polar  coordinate  .system,  and  consists  of  256-1 .4  deg.  radials,  or 
"sectors,”  per  scan.  Due  to  the  requirement  that  the  ASR-9  detect  overhead  microbursts,  it  is  neces¬ 
sary  to  buffer  one  scan’s  worth  of  data  prior  to  the  segment  detection  process  so  that  the  shear  con¬ 
tributions  from  all  azimuths  can  be  taken  into  account.  In  order  to  avoid  having  to  merge  the  shear 
from  opposing  overhead  shear  segments,  which  would  add  an  additional  step  to  the  algorithm,  a 
coordinate  system  shift  is  performed  that  allows  the  algorithm  to  run  on  "composite”  radials  made  up 
of  the  radials  from  0  to  1 80  degrees  and  their  opposing  counterparts.  Because  the  Doppler  signal 
from  opposing  radials  is  of  opposite  sign,  a  coordinate  system  shift  is  needed  to  allow  the  algorithm 
to  run  in  the  same  manner  whether  or  not  it  is  processing  data  overhead. 

Figure  7a  is  an  overhead  view  of  a  microburst  impacting  close  enough  to  the  radar  so  that  some 
of  the  divergent  outflow  passes  overhead.  Figure  7b  plots  the  variation  of  radial  velocity  with  range 
at  the  azimuth  represented  by  the  dotted  line  in  7a.  As  can  be  seen  from  the  solid  line  in  7b,  the  veloc¬ 
ity  signature  contains  a  discontinuity  as  the  range  passes  through  zero.  This  can  be  removed  by  sim¬ 
ply  inverting  the  velocity  on  one  side  of  the  radar,  as  shown  by  the  dotted  line.  The  radial  can  then  be 
processed  in  the  normal  manner,  except  that  the  range  gates  lying  on  the  inverted  side  of  the  compos¬ 
ite  radial  must  be  scanned  in  order  of  decreasing  range.  To  simplify  the  processing,  the  range  of  gate 
values  can  be  thought  of  as  varying  from  -R  to  R  instead  of  the  normal  0  to  R.  The  azimuth  is  varied 
from  0  to  1 80  degrees  instead  of  0  to  360  degrees.  This  is  illustrated  in  Figures  7c  and  7d.  This  tech¬ 
nique  does  not  affect  the  detection  of  non-overhead  events  in  any  way,  and  incurs  only  a  small  pro¬ 
cessing  overhead  (changing  the  sign  of  the  velocity  values  on  one  side  of  the  radar). 

One  additional  issue  that  can  be  of  concern  when  detecting  overhead  events  is  the  quality  of  the 
velocity  data  at  extremely  close  range  (1-2  gates).  At  this  range  the  radar  receivers  can  still  be  recov¬ 
ering  from  the  effects  of  the  transmitted  pulse,  resulting  in  biased  velocity  values.  In  order  to  deter- 
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time:  S2/'0B/06  19:59:54  1  degrees 


Figure  6.  Example  of  a  classicd  microburst  signature.  Velocity  field  is  shown  on  right,  with  corresponding 
r^ectivity  field  on  li^.  The  shear  segments  detected  by  the  microburst  algorithm  are  drawn  in  white. 
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a) 


b) 


Figure  7.  Coordinate  System  Transformation. 


mine  the  magnitude  of  this  problem,  a  study  was  conducted  on  several  overhead  microburst  events 
that  occurred  in  1989  at  the  Kansas  City  field  site.  The  study  found  that  although  the  velocity  values 
for  the  first  two  range  gates  were  sometimes  biased  for  the  lower  reflectivity  events  (<  40  dBZ), 
events  with  higher  reflectivities  had  sufficiently  accurate  velocity  values  to  allow  the  majority  of 
overhead  segments  to  be  detected  in  their  entirety.  Even  in  the  case  where  some  of  the  segments  are 
broken  up  into  two  segments  (one  on  each  side  of  the  radar),  as  long  as  at  least  one  unbroken  segment 
is  detected  the  segment  association  logic  will  still  tend  to  group  all  the  segments  into  the  same  cluster, 
providing  an  increased  tolerance  for  any  short  range  velocity  anomalies.  Based  on  these  findings,  the 
current  version  of  AMDA  has  no  additional  processing  steps  to  correct  for  possibly  biased  velocity 
values  at  short  range.  This  issue  is  still  an  area  of  active  study  however,  as  the  above  findings  are  only 
valid  for  the  transmitter/receiver  chain  they  were  tested  with,  and  the  Lincoln  Laboratory  ASR-9 
testbed  transmitter/receiver  chain  has  been  substantially  modified  since  it  was  moved  from  Kansas 
City  to  Orlando  in  1990, 

Figure  8  illustrates  an  example  of  an  overhead  microburst,  including  the  detected  shear  segments 
and  the  cluster  outline  that  is  displayed  on  the  WSP’s  real-time  displays.  The  velocity  field  indicates 
winds  blowing  away  from  the  radai’  at  all  azimuths,  although  the  velocity  field  is  asymmetrical  be¬ 
cause  the  storm  was  moving  rather  rapidly  to  the  east-northeast. 
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4.3  SEGMENT  DETECTION 


Shear  segment  detection  can  be  broken  down  into  the  following  steps: 

1 .  Detect  the  beginning  of  a  segment, 

2.  Search  along  radial  for  the  end  of  a  segment,  and 

3.  Check  the  segment  for  validity  (length,  shear,  slope). 

4.3.1  DETECT  BEGINNING  OF  SEGMENT 

A  gate  is  flagged  as  possibly  being  the  beginning  of  a  segment  if  the  velocity  at  the  following  gate 
is  higher  than  the  current  gate.  There  is  an  additional  requirement  that  the  increase  between  the  two 
gates  is  not  greater  than  MAX_FIRST_INCREASE  (5.0  m/s).  This  helps  reduce  the  number  of  false 
segments  begun  in  noisy,  low  reflectivity  areas. 

4.3.2  SEARCH  FOR  END  OF  SEGMENT 

The  end  of  a  .segment  is  signalled  by  one  of  two  conditions: 

•  More  than  M  AX_BAD_VALS  ( 2 )  bad  values  encountered.  Bad  values  can 
be  caused  by  clutter  breakthrough,  second-trip  echoes,  and  anomalous  prop¬ 
agation,  and  are  flagged  as  such  by  the  signal  processor  prior  to  AMDA  read¬ 
ing  the  data. 

•  Sum  of  decreases  since  the  most  recent  maximum  value  is  greater  than 
MAX_SUM_DECREASES  (5.0  m/s). 

The  first  condition  is  self-explanatory.  The  second  condition  is  a  little  more  complicated,  and 
represents  the  hcan  of  the  detection  algorithm  since  it  largely  determines  the  overall  length  of  the 
detected  segments.  As  AMDA  is  moving  from  point  to  point  along  the  shear  segment,  a  variable 
sum  jle  creases  is  used  to  monitor  the  decreasing  bend  of  a  velocity  radial.  This  variable  repre¬ 
sents  the  magnitude  of  the  decreasing  bend,  and  will  always  be  a  positive  quantity.  The  following 
rules  are  used  to  increase  and  decrease  the  value  of  the  variable; 

•  Initialize  sum  decreases  to  0.0  at  the  beginning  of  a  segment. 

•  When  a  velocity  decrease  is  found,  the  magnitude  of  the  decrease  is  added 
to  sum  decreases.  When  the  variable  reaches  or  exceeds  MAX_SUM_DE- 
CREASES,  the  growth  phase  of  segment  detection  is  terminated  and  the  seg¬ 
ment  is  b-immed  back  to  the  preceding  maximum  velocity  value. 
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Figure  8.  Detection  of  an  overhead  microburst  event.  Note  that  the  velocity  field  (right)  indicates  winds 
blowing  away  from  the  radar  at  all  azimuths. 
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•  When  a  vckxity  increjse  is  found,  and  sum  decreases  is  non-zero  (recent 
decreasing  trend  encountered),  then  reduce  the  value  of  sum  decreases  by 
SUM_DECREASES_ADJ  ( 1 .0  m/s),  or  set  it  to  0.0  if  the  result  would  be  less 
than  0.0.  The  reason  to  reduce  sum  decreases  by  a  constant  instead  of  the  ac¬ 
tual  magnitude  of  the  velocity  increase  is  to  weigh  the  decreasing  trends  more 
heavily  in  regions  where  the  vel(x:ity  field  is  oscillating  up  and  down.  This 
tends  to  terminate  the  growth  phase  of  noisy  .segments  with  greater  efficien¬ 
cy. 

Figure  9  illustrates  the  segment  growth  process,  including  an  example  of  how  the  algorithm  han¬ 
dles  a  small  anomaly  in  the  velocity  field.  The  ..egment  is  started  at  point  A  because  (  Vb-V^)  >  0.0. 
The  value  oi sumjiecreases  at  point  B  is  iimialized  to  0.0.  It  remains  at  0.0  until  the  first  decrease  is 
detected  (D),  at  which  point  it  is  incrca.scd  by  the  magnitude  of  ( Vc  -  V^).  Although  the  next  point, 
E,  represents  an  increase  with  resnect  to  point  D,  it  is  not  sufficient  to  signal  a  end  to  the  decreasing 
trend  because  the  velocity  is  still  less  than  the  current  maximum  value  (C).  Accordingly,  sum  de- 
creases  is  increased  by  (Vc  -  Vg).  Point  F  represents  a  slight  increase  with  respect  to  point  C,  so 
sum  decreases  is  reduced  by  the  fixed  value  of  1 .0,  to  0.8.  Point  G  is  another  in  crease,  and  the  de¬ 
creasing  tiiend  is  again  reduced  by  1 .0,  wnich  results  in  the  zeroing  of  sumjlecreases  because  it  can¬ 
not  become  negative.  At  this  point,  the  algorithm  has  completely  recovered  from  the  velocity  anom¬ 
aly  at  point  C.  The  segment  continues  to  grow  until  it  reaches  point  I,  which  represents  a  velocity 
decrea.se.  The  value  of  the  decrease  (Vh  -  Vj)  is  added  xo  sumjlecreases,  and  the  process  is  repeated 
for  points  J  and  K,  adding  (Vh  -  Vj)  and  (Vh  -  Vk),  respectively.  At  point  K,  the  segment  growth 
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process  is  terminated  beeausc  j«m_r/eer6’aj^5>5.0ni^s.The  maximum  vekK'ity  value,  point  H  in  this 
ease,  is  taken  as  the  segment  endpoint. 

4.3.3  SEGMENT  THRESHOLDING 

Segment  thresholding  serves  to  remove  those  shear  .segments  that  arc  t(X)  short  or  ttx)  weak  to 
have  been  caused  by  a  significant  divergent  outflow.  The  .segment  thresholding  steps  are: 

•  Segment  length  must  be  greater  than  or  equal  to  MIN_SEG_LENGTH  (().a6 
km). 

•  Segments  must  exhibit  a  differential  velocity  greater  than  or  equal  to 
M1N_SEG_SHEAR  (7.8  m/s). 

•  The  shear  slope  at  each  point  in  the  segment  mustcxcced  MIN_SEG_SLOPE 
(1.8  m/s  /  km).  A  dual-window  approach  utilizing  both  an  1 1-point  and  a 
7-point  window  (centered  at  the  test  point)  is  used  to  compute  the  slope,  ex¬ 
cept  near  the  endpoints  of  the  segment  where  a  shorter  single  window  (9, 7, 

5,  or  3-point)  is  used  due  to  lack  of  sufficient  data.  The  slope  for  each  of  the 
two  windows  is  computed  using  the  velocity  difference  between  the  window 
endpoints,  and  the  larger  of  the  two  results  is  taken  as  the  final  slope  value. 

The  two  windows  aie  necessary  to  handle  the  cases  of  longer  segments  wnere 
small  dips  or  ri.ses  in  a  region  of  relatively  low  shear  can  cause  a  single  win¬ 
dow  scheme  to  falsely  reject  a  segment.  This  is  generally  not  the  case  near 
the  segment  endpoints,  and  the  use  of  only  a  single  window  in  the.se  regions 
has  not  had  a  negative  impact. 

If  the  slope  at  any  point  is  less  than  the  threshold,  the  segment  is  split  at  that 
point.  Any  segment  that  is  split  up  into  two  or  more  segments  must  pas.> 
through  the  minimum  length  and  shear  tests  once  more  to  ensure  segment  va¬ 
lidity. 

The  cuiTent  value  of  1.8  m/s/  km  for  MIN_SEG_SLOPE  is  somewhat  lower 
than  the  2.5  m/s/  km  slope  as  calculated  from  the  definition  of  a  Microburst 
(10  m/s  over  4  km).  A  setting  lower  than  2.5  is  necessai^  to  handle  regions 
of  shallow  slope  such  as  illustrated  in  Figure  2b,  and  the  value  of  1 .8  has  been 
experimentally  determined  to  reject  excessively  large  regions  of  shallow 
slope  while  retaining  the  regions  that  naturally  occur  within  a  microburst. 

The  slope  check  thresholding  test  is  the  most  computationally  intensive  of  the  three  thresholding 
steps,  and  for  that  reason  the  segment  length  and  strength  tests  are  mn  first.  The  great  majority  of 
invalid  segments  will  be  filtered  out  by  the  first  two  tests,  thereby  reducing  the  number  of  segments 
that  requiie  the  slope  check. 
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4.4  SEGMENT  TEMPORAL  SMOOTHING 

The  ASR-9  antenna  seans  onec  every  4.8  sceonds,  providing  an  opportunity  to  temporally 
sm(X)th  the  shear  segments  to  decrease  the  scan-to-scan  variability.  In  the  AMDA  implementation, 
evci7  .segment  has  ass(x:iatcd  with  it  a  variable  called  seg_a}>e,  which  stans  out  at  a  value  of  1  when 
the  segment  is  first  detected.  For  each  successive  .scan  on  which  the  segment  is  present  (determined 
by  overlap),  segjjge  is  incremented  by  I,  until  it  reaches  a  maximum  value  of  iviAX_SEG_AGE 
( 7  ).  For  each  .scan  where  the  segment  is  missing,  seg_age  is  decreased  by  1 .  A  segment  is  automati¬ 
cally  in.scrtcd  (coasted)  into  a. scan  it'seg_age>=  MIN_SEG_AGE_INSERT(  3  )and  the  segment  is 
missing  from  the  cunent  .scan. 

Figure  10  shows  the  results  from  a  Monte-Carlo  simulation  of  the  above  algorithm  for  three  dif¬ 
ferent  "segment  present”  percentages.  From  the  graph,  it  can  be  seen  that  for  segments  present  on 
only  25  percent  of  the  scans,  the  probability  that  seg  age  reaches  3  is  quite  small,  resulting  in  the 
in.scrtion  of  coast  segments  roughly  3  percent  of  the  time.  For  segments  that  are  present  50  percent  of 
the  time,  coast  segments  will  be  created  58  percent  of  the  time,  and  when  segments  exist  on  75  per¬ 
cent  of  the  scans,  a  coast  segment  will  be  inserted  98  percent  of  the  time.  The  segment  age  require¬ 
ment  of  3  was  chosen  because  of  the  combination  of  a  small  percentage  of  ‘false’  in.sertions,  and  a 
high  percentage  of  "true’  inseilions. 


X  Segment  present  75%  of  time 
A  Segment  present  50%  of  time 
□  Segment  present  25%  of  time 


Figure  10.  Probability  Distribution  of  Parameter  segage. 
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The  above  argument  is  directed  toward  the  case  where  a  given  segment  is  flickering  on  and  off 
because  it  is  near  one  of  the  thresholds.  Another  ca.se  to  consider  is  that  of  a  point  target  (aircraft) 
causing  a  segment  miss  on  several  consecutive  scans.  Assuming  that  the  segment  had  been  fairly 
steady  just  prior  to  the  point  target  andicjEja^'c  had  reached  the  maximum  value  of  7,  a  coast  segment 
would  be  created  for  up  to  four  (7-3)  scans  in  which  the  segment  was  missing.  This  conesponds  to 
4  X  4.8  seconds/scan  =  19.2  seconds  of  enhanced  segment  continuity,  which  is  adequate  when 
compared  with  the  time  a  typical  moving  point  tiuget  will  remain  within  a  radial. 

4.5  SEGMENT  ASSOCIATION 

Segment  asscx'iation  is  done  in  a  manner  similar  to  the  TDWR  algorithm.  Segments  on  neighbor¬ 
ing  radials  arc  tested  for  range  overlap,  and  grouped  together  into  "clusters”  when  overlap  occurs. 
There  are  three  main  steps  in  the  grouping  process. 

1 .  Using  an  "azimuthal  association  window,”  determine  if  two  segments  are 
close  enough  in  azimuth  to  be  subjected  to  the  range  overlap  test. 

2.  If  so,  begin  a  new  cluster  when  the  segments  also  overlap  each  other  in  range. 

3.  Continue  to  grow  the  cluster  until  no  more  overlapping  segments  are  encoun¬ 
tered  within  the  azimuthal  association  window. 

The  TDWR  algorithm  utilizes  a  fixed  azimuthal  association  window  (usually  2  deg.)  when  per¬ 
forming  segment  as.s(K'iation.  In  the  case  of  the  ASR-9,  which  must  detect  microbursts  at  close  range 
owing  to  on-airpoii  location,  nearby  point  targets  or  ground  clutter  can  often  obscure  a  larger  azi¬ 
muthal  region  than  tho.se  at  greater  distances.  This  situation  is  depicted  in  Figure  1 1  a.  To  address  this 
problem,  the  AMDA  approach  is  to  make  the  azimuthal  association  window  a  linear  function  of 
range  (The  range  is  determined  from  the  midpoint  of  the  segment  under  consideration ).  The  maxi¬ 
mum  and  minimum  limits  of  the  window  size  are  specified  as  variable  site  pai  ameters  (See  Appen¬ 
dix  A  for  full  parameter  names).  The  maximum  window  size  is  typically  set  at  15.4  deg.  at  zero 
range,  ramping  down  to  a  minimum  value  of  4.2  deg.  at  a  range  of  10.0  km  (Figure  1 1  b).  These  pa¬ 
rameters  are  generous  when  compared  with  the  TDWR  setting  of  2  deg.,  and  it  is  true  that  these  set¬ 
tings  will  result  in  some  clusters  containing  only  widely  scattered  segments  at  close  ranges.  There  is 
an  additional  AMDA  processing  step,  segment  density  thresholding,  to  address  this  issue,  and  the 
explanation  is  defeiTed  to  that  section. 

The  overlap  test  simply  tests  if  the  inner  SEG_OVERLAP_PERCENTAGE(  50  %)  of  neighbor¬ 
ing  segments  have  any  overlap.  During  AMDA  optimization  experiments,  this  overlap  requirement 
has  been  varied  between  25  percent,  which  tended  to  fragment  segments  from  the  same  microburst 
into  different  clusters,  and  75  percent,  which  tended  to  group  segments  from  multiple  microburets 
into  a  single  cluster.  The  value  of  50  percent  appears  to  be  a  good  compromise  between  these  two 
extremes. 

To  be  considered  valid,  clusters  must  contain  at  least  MIN_CLUSTER_SEGS  ( 4 ).  If  this  is  true, 
then  the  cluster  is  passed  on  to  the  reflectivity  thresholding  module. 


a) 


b) 


Figure  II.  Association  Criteria  Range  Dependence. 


4.6  REFLECTIVITY  THRESHOLDING 

The  ASR-9  dual-beam  velocity  estimator  may  not  reliably  predict  the  low-level  radial  velocity 
field  in  all  meteorological  situations.  Analysis  of  data  acquired  in  Huntsville,  AL  and  Kansas  City, 
KS  revealed  that  shear  segments  located  primarily  on  the  fringes  of  storm  cells  were  a  major  source 
of  false  alarms,  most  often  due  to  errors  in  the  dual-beam  velocity  field  at  these  locations.  In  the 
majority  of  cases,  overhanging  precipitation  had  contaminated  the  power  spectrum  in  excess  of  the 
dual-beam  algorithm’s  ability  to  compensate, causing  the  erroneous  velocity  values.  An  example  of 
this  type  is  presented  in  the  performance  evaluation  section  of  this  report.  Since  the  situation  occurs 
most  frequently  at  the  edges  of  storm  cells,  AMDA  includes  a  reflectivity-based  thresholding  step  to 
filter  out  this  type  of  false  alarm. 

The  reflectivity  thresholding  process  is  illustrated  in  Figure  1 2.  First,  segments  with  greater  than 
REFL_THRESH_LENGTH_PERCENT  (50%)  of  their  length  lying  within  “core”  reflectivity 
(  >  REFL_THRESH_CORE_DBZ  (40  dBZ) )  are  marked  as  valid  microburst  segments.  Next,  to 
reduce  the  chances  of  discarding  shear  segments  associated  with  a  true  microburst,  any  other  seg¬ 
ments  that  lie  within  REFL_THRESH_DIST  (0.5  km)  of  a  valid  segment,  as  measured  between  the 
two  segment  mid-points,  are  also  marked  as  valid.  All  other  segments  in  the  cluster  are  discarded. 

The  reflectivity  thresholding  process  significantly  reduces  the  false-alarm  rate  in  the  environ¬ 
ments  studied  thus  far.  A  typical  case  illustrating  the  effect  of  the  thresholding  is  shown  in  Figure  1 3. 
In  this  case,  a  divergent  signature  is  clearly  present  in  the  velocity  field  (right)  and  a  cluster  of  seg- 
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Core  Reflectivity 


Valid  Shear 
Segments 


Fringe 

Reflectivity 


Invalid  Shear  Segments 
(Less  than  50%  of  segment  length 
lies  within  core  reflectivity) 


Figure  12.  Kefleaivity  Thresholding  oj  Shear  Segments. 


mcnts  (6)  that  passes  the  area  and  str  ength  thresholds  is  present.  However,  examination  of  the  reflec¬ 
tivity  field  (left),  indicates  that  the  majority  of  the  segments  are  located  in  a  region  of  low  reflectivity, 
and  are  potentially  invalid.  Four  out  of  the  six  segments  were  in  fact  rejected  by  the  reflectivity 
threshold,  and  the  remaining  segments  were  insufficient  to  qualify  as  a  microburst.  Subsequent  anal¬ 
ysis  of  this  case  showed  the  receding  portion  of  the  velocity  signature  was  being  caused  by  precipita¬ 
tion  aloft,  and  the  reflectivity  threshold  had  successfully  eliminated  a  false  alarm. 

In  a  dry  environment,  the  outcome  of  using  this  method  remains  to  be  seen,  and  the  algorithm 
will  probably  need  modification  in  order  to  avoid  the  rejection  of  dry  events  occurring  simultaneous¬ 
ly  with  events  containing  higher  reflectivities.  One  possibility  currently  being  evaluated  is  to  deter¬ 
mine  the  core  reflectivity  threshold  on  a  storm-cell  by  storm-cell  basis  instead  of  utilizing  a  single 
threshold  for  the  entire  scan.  Improvements  to  the  dual-beam  algorithm  are  also  expected  to  reduce 
the  overall  need  for  reflectivity  thresholding. 

4.7  SEGMENT  DENSITY  THRESHOLDING 

In  the  description  of  the  initial  segment  association  step,  we  mentioned  that  the  relatively  large 
radial  gap  allowed  between  associated  segments  could  result  in  large  clusters  with  only  a  sparse 
number  of  segments.  There  are  actually  four  cases  to  consider. 

The  first  case,  shown  in  Figure  14a,  is  that  of  the  ideal  cluster,  with  no  missing  segments.  This 
case  needs  to  be  passed  unaltered  through  the  segment  density  thresholding  process.  The  second 
case.  Figure  14b,  is  that  of  a  cluster  with  multiple  missing  segments,  but  with  a  high  density  of  con¬ 
tiguous  detections  in  the  regions  neighboring  the  missing  segment  gap.  The  most  common  causes  of 
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Figure  13.  R^ctivity  thresholding  succes^ully  elumnating  a  false  alarm  on  the  fringe  of  a  storm  cell.  The 
group  of  segments  in  this  case  met  all  the  microburst  area  strength  requirements,  but  failed  the  reflectivity 
threshold. 
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a)  Ideal  cluster  with  no 
missing  segments 


b)  Cluster  with  multiple  missing  segments, 
but  high  neighboring  segment  density. 


C)  Cluster  with  valid  high-density  region 
and  invalid  low-density  region; 
reject  low-density  portion 


d)  Cluster  with  sparse  segment  density; 
reject  entire  duster 


Figure  14.  Typical  Segment  Density  Distributions  and  Associated  Clustering  Decisions. 


this  condition  arc  ground  clutter  and  point  target  contamination.  It  is  desirable  to  accept  cases  such  as 
this  as  a  valid  cluster.  The  third  case.  Figure  1 4c,  consists  of  a  cluster  with  a  region  of  high  segment 
density,  and  also  a  low  density  region.  The  desired  action  in  this  case  is  to  trim  off  the  low  density 
region  while  retaining  the  high  density  region.  The  fourth  case.  Figure  14d,  is  that  of  a  cluster  that 
contains  a  very  sparse  disti  ibution  of  segments,  in  which  case  the  entire  cluster  should  be  rejected. 

To  achieve  the  desired  density  thresholding  behavior,  AMDA  utilizes  a  sliding  window  of  width 
[(AZIMUTH_G  AP(Range)  x  2)  +  1  ]  to  compute  the  segment  density  at  each  radial.  The  azimuth  gap 
is  the  same  as  that  used  in  the  segment  association  process,  and  is  a  function  of  range.  The  sliding 
window  is  always  centered  on  the  segment  being  evaluated,  so  the  window  size  near  the  cluster  edges 
must  be  reduced  if  necessary  to  avoid  extending  beyond  the  cluster  azimuth  limits.  For  example,  the 
window  size  for  the  second  and  third  azimuth  positions  in  a  cluster  is  limited  to  three  and  five  radials, 
respectively.  No  test  is  performed  on  the  edge  segments  themselves.  The  segment  density  is  simply 
calculated  as  the  percentage  of  segments  present  in  the  window,  and  must  be  greater  than  or  equal  to 
SEG_DENSITY_THRESH  (50%)  or  the  cluster  is  split  at  that  point.  All  resulting  sub-clusters  are 
then  required  to  meet  the  minimum  segment  requirement  of  CLU_MIN_SEGS(4),  or  they  are  dis¬ 
carded. 

4.8  CLUSTER  SHEAR  DETERMINATION 

Even  with  the  spatial  filtering  and  data  consistency  checks  performed  up  to  this  point,  it  is  still 
possible  to  encounter  occasional  shear  segments  with  spurious  values.  This  is  obviously  undesirable 


when  estimating  the  overall  shear  c^f  the  elustcr,  because  the  overall  shear  is  nomially  defined  as  the 
maximum  value  aeross  all  shear  segments,  and  a  single  outlier  can  dramatically  influence  the  maxi¬ 
mum  value.  On  the  other  hand,  snKX>thing  the  shear  values  excessively  will  result  in  the  loss  of  the 
true  detected  maximum  shcai.  AMDA  uses  a  3-point  median  test  tocheck  through  the  segments,  and 
replace  any  detected  outliers  with  the  median  of  the  shear  from  the  sunounding  three-point  window 
(Figure  15).  This  process  will  remove  most  outliers,  while  retaining  the  true  maximum  shear  value  in 
the  majority  of  cases.  Subsequent  to  the  median  filter,  the  cluster  sheai’  is  set  equal  to  the  maximum 
sheai’  value  acro.ss  all  segments  in  the  cluster. 


If  shear  value  of  a  segment  is  more  than  20  percent  higher  than  the  median 
of  the  three-segment  window,  replace  shear  value  with  median  value. 


Figure  15.  Shear  Magnitude  Outlier  Rejection. 


4,9  AREA  THRESHOLDING 

The  purpose  of  area  thresholding  is  to  remove  clusters  of  .segments  which,  because  of  their  small 
size,  are  probably  not  true  microburst  events.  This  is  not  a  trivial  problem,  as  some  of  the  stronger 
microbursts  are  initially  vei7  compact,  and  it  is  essential  to  avoid  filtering  out  these  microbursts 
while  attempting  to  improve  the  overall  false-alann  rate.  To  help  alleviate  the  problem,  it  is  useful  to 
note  the  relationship  between  Probability  of  False  Alarm  (PEA)  and  shear  magnitude.  Figure  16  is  a 
plot  of  PFA  vs.  shear  magnitude  using  the  average  of  the  1991  and  1992  scoring  results  presented  in 
Section  5  of  this  report.  As  seen  in  the  figure,  the  likelihood  of  a  cluster  being  a  false  alarm  decreases 
rapidly  with  increasing  shear  magnitude.  This  information  can  be  used  to  advantage  in  the  area 
thresholding  process  by  using  a  threshold  that  is  inversely  proportional  to  shear  magnitude. 

A  second  consideration  is  that  of  asymmetrical  microbursts,  when  the  percentage  of  microburst 
area  detected  can  sometimes  decrease  with  decreasing  range.  As  an  asymmetrical  microburst  (with 
major  axis  perpendicular  to  the  radar  beam)  nears  the  radar,  the  viewing  angle  of  the  radar  renders 
some  of  the  velocity  signature  of  the  microburst  invisible,  due  to  the  geometry.  This  is  illustrated  in 
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Figure  16.  ProbabUiry  of  False  Alarm  rv.  Shear  Magnitude. 


Figure  1 7.  When  the  microburst  is  relatively  distant  (left),  the  radar  “secs”  almost  the  entire  region  of 
shear,  but  as  the  microburst  moves  closer  to  the  radar  (right),  the  radial  component  of  the  shear  de¬ 
creases  at  the  wider  angles,  and  some  of  the  segments  are  not  detected  because  the  shear  falls  below 
the  detection  threshold. 

The  previous  considerations  can  be  utilized  effeetively  by  making  the  area  threshold  a  function 
of  both  cluster  shear  strength  and  range.  The  relationship  used  by  AMD  A  is  shown  in  Equations  1  to 
3.  The  parameters  used  in  the  equations  are  defined  as  follows  (current  settings  in  parentheses); 

Areaniin  Absolute  minimum  area  threshold  (0.5  km^) 

Rniin,  Rniax  Minimum  and  maximum  range  boundaries  (  2.0, 16.0  km  ) 

Shearniin,  Shearmax  Minimum  and  maximum  shear  boundaries  (  7.8,  27.5  m/s  ) 

Rsiope  Rate  of  change  of  area  threshold  with  respect  to 

range  (O.lkm^/km) 

A,  B,  C  Area  threshold  shear  correction  coefficients  (30.0,  2.5,  -1.0) 


Portioiu  of  asymmetrical  or 
li"  microbursts  can 
become  less  visible  to  the 
radar  as  the  range 
decreases,  due  to  viewing 
angle. 


Figure  17.  Detected  Microburst  Area  vs.  Range. 
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Note  that  the  parameter  names  are  not  the  full  symbolic  names  (i.e.  CLU_M1N_AREA)  used  in  the 
algorithm  parameter  tile.  The  eomplete  list  of  symbolic  parameter  names  can  be  found  in  Appendix 
A. 

The  area  threshold  for  a  given  cluster  of  segments  is  computed  using  the  sum  of  an  absolute 
minimum  area  threshold  and  two  correction  terms; 

=  Arc'a„,„,  +  AAreara.,se  +  ‘AArea^,^,,,  (1) 

The  range  coirection  term  is  computed  using  a  simple  linear  relationship  with  one  adjustable 
parameter,  Rsiopc  iind  has  the  effect  of  increasing  the  area  threshold  vs.  range. 

AAreUran^e  ~  ^slope^^  ~~  [2] 

The  shear  correction  term  is  calculated  using  a  non-linear  equation  designed  to  match  the 
shape  of  the  PFA  vs.  Shear  Magnitude  cui've.  This  temi  is  small  for  strong  events,  but 
increases  rapidly  as  the  event  strength  decreases. 

131 

The  two  collection  terms  are  both  prevented  from  becoming  negative  or  excessively  lai'ge  by 
bounding  the  independent  variable  in  each  case.  The  limits  are  shown  below. 

R,  Rmi,.  <  R  <  Rmax  Sheai'min  <  Shcai'  <  Shearnvw 

Rmin.  R  <  Rtnin  Shear  =  ,  Shearmi,,.  Sheai-  <  Shearmi,, 

Rmaxi  R  ^  Rrmx  Shcaimax'  Shcai  >  Shcainux 

k 

The  parameter  settings  shown  above  were  determined  in  an  iterative  manner,  whereby  single 
paiameters  were  vai  ied  in  turn,  and  the  algorithm  scoring  results  were  examined.  At  a  range  of  2km, 
the  current  settings  result  in  an  area  threshold  that  ranges  from  0.5  km^  for  very  strong  events  (>27.5 
m/s)  to  2.5  km^  for  the  weakest  events  (7.5  m/s).  A  contour  plot  illustrating  the  operational  area 
thresholds  (boldface  -  in  square  kilometers)  is  provided  in  Figure  18.  Because  the  microburst  area 
distribution  is  sensitive  to  changes  in  other  parameters,  especially  the  allowed  radial  gap,  the  area 
thresholding  parameters  must  be  re-optimized  when  other  paiameters  aie  changed. 

4.10  CLUSTER  TEMPORAL  SMOOTHING 

To  provide  a  more  consistent  microburst  appearance  from  scan  to  scan,  a  simple  temporal  filter  is 
utilized.  To  prevent  spurious  clusters  from  being  output  as  valid  microbursts,  a  cluster  must  be  pres¬ 
ent  on  three  successive  scans  to  be  declared  valid.  This  introduces  a  10  second  latency  period  in  the 
initial  detection  phase,  but  because  microbursts  are  almost  always  detected  in  their  growing  stage, 
this  amount  of  latency  is  not  viewed  as  a  serious  problem.  To  reduce  the  amount  of  microburst  “flick¬ 
ering”  that  can  occur  while  an  event  is  dying  out,  microbursts  are  “coasted”  for  a  maximum  of  12 
scans  ( 1  minute).  This  tends  to  increase  the  overall  false-alarm  rate  slightly  because  false  alarms  are. 
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EVENT  STRENGTH  (M/S) 


Figure  18.  AM  DA  Area  Thresholding. 


in  effect,  “amplified,”  but  this  minor  FAR  increase  is  preferable  to  providing  ATC  with  an  inconsis¬ 
tent  display. 


-31- 


4. 1 1  SHAPE  ALGORITHM  AND  FINAL  OUTPUT 


The  temporal  smoothing  of  clusters  is  the  final  step  of  the  microburst  algorithm  that  is  specific  to 
the  ASR-9.  From  that  point,  microburst  warnings,  along  with  the  segment  clusters  for  each  micro¬ 
burst,  are  sent  to  the  shape  algorithm  [Wilson,  1992]  running  on  a  Sun  Workstation.  This  algorithm, 
identical  to  the  algorithm  used  by  TDWR,  creates  “band-aid”  shape(s)  for  each  microburst,  which 
arc  then  sent  with  the  microbui-st  shear  value  to  the  Geographical  Situation  Display  (GSD)  in  the 
Control  Tower.  An  example  of  the  GSD  display  is  .shown  in  Figure  19. 
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Figure  19.  Geographic  Situation  Display.  Filled  red  circles  are  detected  microbursts  with  the  measured 
wind  changes  across  the  event  (in  knots)  indicated.  Open  red  circles  show  divergent  outflows  with  mea¬ 
sured  wind  changes  less  than  30  knots  ( "wind  shear  with  loss”).  Purple  arc  is  a  detected  gust  front  with 
the  estimated  location  of  the  front  10  and  20  minutes  in  the  future  indicated  by  dashed  purple  lines;  an 
estimate  of  the  wind  speed  and  direction  behind  the  front  is  given  by  the  purple  vector.  Six  levels  of pre¬ 
cipitation  r^ectivity  are  shown  and  the  speed  and  direction  of  storm  movement  is  given  by  the  black 
vector  and  associated  number.  The  blue  arrows  and  numbers  show  the  wind  speed  and  direction  mea¬ 
sured  by  LIMAS  anemometers. 
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5.  MICROBURST  ALGORITHM  PERFORMANCE  EVALUATION 


5.1  EVALUATION  STRATEGY  AND  RESULTS 

Signal-processing  and  meteorological  detection  algorithms  for  the  ASR-9  WSP  have  been  de¬ 
veloped  and  validated  using  data  collected  in  Huntsville,  AL.  Kansas  City,  KS  and  Orlando,  FL  with 
the  Lincoln  Laboratory  testbed.  This  section  discusses  the  pciformance  of  the  “mature”  version  of 
the  ASR-9  Microburst  Detection  Algorithm  (AMDA),  using  data  from  eleven  active  weather  days 
during  the  1991  Orlando  operational  test,  and  ten  active  weathei  days  during  the  1992  test.  AMDA 
output  was  scored  against  corresponding  "truth”  data  from  the  TDWR  pencil-beam  radar  located 
approximately  eight  kilometers  south  of  the  Orlando  ASR  testbed.  Truth  data  were  generated  by  ex¬ 
amining  TDWR  low-level  PPI  scans  taken  at  one-minute  intci-vals  during  the  active  periods.  Using 
truthing  software,  locations  (polygonal  outlines)  and  intensities  of  microbursts  were  entered  into  a 
computer  file.  An  automated  scoring  program  then  compared  the  locations  of  all  ASR  alarms  within 
30  seconds  of  a  truth  time  with  the  locations  of  tme  events  at  that  time,  scoring  an  alamt  a  "hit”  if  it 
overlapped  a  truth  region,  and  a  "false  alarm”  if  it  did  not.  Truth  events  that  were  not  intersected  by 
an  ASR  alarm  were  scored  as  "misses.”  Since  the  majority  of  TDWR  low-level  scans  were  "sector” 
scans  centered  over  the  Orlando  International  Airpon,  azimuthal  limits  were  included  in  the  truth 
files  and  applied  by  the  automated  scoring  program  so  that  only  aUums  within  the  region  of  available 
ti\ith  data  were  scored.  An  alarm  which  lay  panially  within  the  azimuthal  bounds  was  not  scored 
unless  it  intersected  a  truth  outline,  in  which  ca.se  it  was  considered  a  hit. 

Scoring  results  for  1991  and  1992  ai  e  presented  in  Tables  1  and  2.  It  is  important  to  note  that  this 
end-to-end  evaluation  ofWSP  microburst  detection  capability  was  highly  dependent  upon  the  accu¬ 
racy  of  base  data  generated  by  the  signal  processing  algorithm.  Probability  of  Detection  (POD)  and 
Probability  of  False  Alarm  (PFA)  are  Indicated  as  a  function  of  microburst  intensity  (DELTAV)  in 
meters/second.  The  statistics  are  cumulative  from  right  to  left.  Denominators  in  the  two  ratios  were, 
respectively,  the  total  number  of  truth  observations  and  ASR  alaims  in  the  scored  data  set.  Thus,  for 
microbursts  located  within  the  operationally  significant  near-airport  zone  extending  8  kilometers  (5 
miles)  from  the  ASR  and  possessing  maximum  deltav  of  at  least  1 .5  m/s  (30  knots),  the  ASR-WSP 
system  was  greater  than  95  percent  successful  in  detecting  the  events.  The  corresponding  PFA’s  were 
seven  and  1 0  percent,  respectively,  for  the  two  years.  Significantly  larger  false-alarm  probabilities 
occur  for  divergent  outflows  weaker  than  1 5  m/s.  Fonunately,  false  alarms  in  this  category  are  of  less 
concern  than  those  for  stronger  outflows  because  1 )  most  pilots  will  choose  to  continue  a  landing  or 
take-off  when  given  a  “wind  shear  with  loss”  alert  generated  by  these  weaker  events,  and  2)  Orlando 
microbursts  are  often  slightly  asymmetric,  causing  legitimate  discrepancies  in  observed  velocity  dif¬ 
ferentials  between  the  ASR  and  truth  radar.  In  the  latter  case,  false  alarms  and  microburst  misses 
occur  frequently  if  one  radar  observes  weak  microburst-strength  shear  (slightly  above  1 0  m/s  within 
4  km)  while  the  other  radar  observes  below-microburst-thrcshold  shear.  Another  statistic  provided 
in  Tables  1  and  2,  Shear  Ratio  (SR),  was  the  mean  difference  in  Intensity  of  the  true  events  as  seen  by 
the  ASR-WSP  system  versus  the  truth  radar.  Thus,  for  microbursts  stronger  than  25  m/s  (50  knots), 
the  ASR-WSP  system  typically  underestimated  the  maximum  divergence  by  approximately  5  m/s 
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( 10  knots).  It  should  be  noted,  however,  the  AMDA  measure  ot  maximum  divergenee  (X'eurs  along  a 
single  azimuth,  but  the  human  "tr'-thcr”  might  consider  multiple  azimuths  when  assigning  strengths 
to  truth  events. 


Table  1. 

Hit-Miss  Scoring  Results  for  1991  Orlando  Test 


1  PERFORMANCE  STATISTICS  BY  DELTAV  (m/s)  OU  T  TO  8  km 

>=  10 

>=  15 

>=  20 

>=  25 

POD 

5447/  5900=0.92 

4026/4129=0.98 

2424/2  429=1.00 

1051/  1052=1.00 

FAR 

1343/  8586=0.16 

427/  5863=0.07 

46/  2865=0.02 

0/  1040=0.00 

SR 

1.44 

-0.40 

-2.23 

-4.11 

Table  2. 

Hit-Miss  Scoring  Results  for  1992  Orlando  Test 


1  PERFORMANCE  STATISTICS  BY  DELTAV  (m/s)  OUT  TO  8  km 

>=  10 

>=  15 

>=  20 

>=  25 

POD 

531 7  /  6079=0.87 

3776/  3935=0.96 

2018/  2067=0.98 

1011/  1  037=0.97 

FAR 

2251/  8534=0.26 

489/  4981=0.10 

5/1  915=0.00 

0/  617=0.00 

SR 

0.78 

-0.82 

-3.33 

-5.46 

An  alternate  .scoring  technique,  known  as  path-oriented  scoring,  utilized  dual-Doppler  wind 
measurements  to  indicate  the  wind  shear  loss/gain  that  would  be  experienced  by  aircraft  along  the 
runways  and  their  respective  approach  or  departure  corridors.  The  TDWR  and  University  of  North 
Dakota  (UND)  radars  were  the  sources  of  the  dual-Doppler  truth.  This  type  of  scoring  provided  a 
more  critical  evaluation  of  ASR-WSP  performance  because  favorable  scoring  results  required  high¬ 
ly  accurate  detection  of  event  location  and  strength  along  these  nan  ow  "corridors.”  The  losses  indi¬ 
cated  by  the  truth  were  compared  to  those  issued  by  the  ASR-WSP  system  on  nine  days  during  1991 
operations,  and  yielded  the  results  in  Table  3.  Here,  PCL^ssILoss)  represents  the  probability  the 
ASR-WSP  issued  a  loss  alert  (wind  shear  with  loss  less  than  30  knots)  or  microburst  warning  (mi¬ 
croburst  with  loss  greater  than  30  knots)  for  a  true  loss  of  less  than  30  knots.  Similarly,  P(LosslMB) 
indicates  the  probability  the  ASR-WSP  issued  a  loss  alert  or  microbursi  warning  given  there  existed 
a  loss  greater  than  30  knots,  and  P(MBIMB)  gives  the  probability  a  mic’’oburst  warning  was  declared 
for  these  same  conditions.  PFA(MB)  was  the  probability  the  ASR-WSP  issued  a  microburst  warning 
when  no  loss  was  evident  in  the  dual-Doppler  truth.  PFA(Loss)  gives  the  probability  the  ASR-WSP 
declared  a  microburst  warning  or  loss  alert  under  the  same  circumstances.  P(OW),  or  "overwam- 
ing,”  was  the  probability  a  microburst  warning  was  given  when  dual-Doppler  data  revealed  a  loss 
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alert  was  more  appiopriatc.  These  lesulls  are  similar  to  those  produeed  by  the  hit-miss  stiategy. 
Slightly  lower  deteetion  probabilities  resulted  from  the  aforementioned  ASR-WSP  underestimates 
of  the  stronger  diveigenees,  and  by  the  more  rigorous  requirements  of  this  scoring  method.  Differ¬ 
ences  in  the  twc)  scoring  method  statistics  emphasize  leniencies  in  the  hit-miss  scoring  mcth(xl,such 
as  rewarding  full  detection  for  partially-observed  or  grossly  over  or  underestimated  (intensity) 
events. 


Table  3. 

Path-Oriented  Scoring  Results  for  1991  Test 


Detection  Probabilities 

False-Alarm  Probabilities 

P  (Loss/Loss) 

0.72 

PFA  (MB) 

0.05 

P  (Loss/MB) 

0.97 

PFA  (Loss) 

0.10 

P  (MB/MB) 

0.84 

P(OW) 

0.26 

5.2  INVESTIGATION  OF  ALGORITHM  FAILURES 

Oecunenccs  of  ASR-WSP  false  alarms  and  missed  mierobursts  in  the  scored  data  were  identi¬ 
fied  by  overlaying  single-Doppler-based  truth  polygons  and  ASR-WSP  alarms  over  images  of  the 
coiresponding  ASR  and  TDWR  (or  UND)  surface  radar  data.  This  practice  often  led  to  discovery  of 
the  cause  of  a  false  alarm  or  miss.  For  example,  an  ASR-WSP  false  alarm  bounding  a  region  of  di¬ 
vergence  clearly  visible  in  the  .ASR  velocity  field,  but  not  evident  in  the  corresponding  region  of  the 
nearest  TDWR  velocity  data,  prompted  a  comparison  of  the  ASR  and  TDWR  retlectivity  fields 
within  the  alai  m  region,  often  revealing  substantially  larger  areal  coverage  of  precipitation  in  the 
ASR  data.  This  suggested  the  five-degree- wide  beam  of  the  ASR  was  observing  elevated  reflectiv¬ 
ity  structures  in  these  cases,  and  this  was  verified  by  subsequently  examining  true  or  synthesized 
RHl  data  through  the  area  of  interest.  In  the  color  plots  that  follow,  alarms  generated  by  the  ASR- 
WSP  are  represented  by  white  polygons,  and  divergence  segments  (azimuthal  lines  along  which  sig¬ 
nificant  divergence  existed  in  ASR  data)  ai'e  represented  by  white  segments.  Locations  of  trath 
events  are  indicated  by  red  polygons.  In  Figure  20,  a  fast-moving  gust  front  (note  velocity  folding  in 
C-band  TDWR  data)  spawned  by  storms  west  of  the  airport  initiated  new  cell  growth  southeast  of 
the  ASR,  immediately  behind  the  frontal  boundary.  Velocity  data  show  the  ASR  (upper  right)  ob¬ 
served  a  22.7  m/s  maximum  loss  centered  seven  kilometers  east-southeast  of  the  radar,  while  the 
TDWR  (lower  right)  saw  only  slight  divergence  at  that  location.  Of  particular  note  is  the  presence  of 
approaching  (negative)  velocities  within  the  alarm  region  in  the  ASR  data.  A  comparison  of  the  re¬ 
flectivity  fields  (left-hand  panels)  associated  with  the  alarm  shows  far  less  signal  return  in  the 
TDWR  data  (lower  left).  A  synthesized  RHI  scan  (Fig.  21)  created  from  TDWR  volume  scan  data 
through  the  alarm  region  revealed  easterly  winds  existed  within  an  elevated  reflectivity  core  (si¬ 
tuated  above  folded  velocities  associated  with  the  fast-moving  gust  front),  while  surface  winds  were 
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westerly  clue  to  the  gust  frontal  passage.  The  ASR  was  sensitive  to  the  easterly  (inbound)  motion 
becau.se  of  its  wide  beamwidth.  Thus,  this  is  an  example  in  which  physical  characteristics  of  the  ASR 
led  to  issuance  of  false  microburst  warnings  by  the  microburst  detection  algorithm. 


The  ASR-9  Microburst  Detection  Algorithm’s  use  of  "noisy”  vel(xity  estimates  in  regions  of 
low  signal  return  also  led  to  issuance  of  false  alarms.  Large-magnitude  velocity  estimates  some¬ 
times  (x;cuiTed  in  ASR  range  bins  in  areas  of  low  signal  rciui-n,  and  the  incorporation  of  these  esti¬ 
mates  by  the  detection  algorithm  prcxluced  fictitiously  high  divergences.  Figure  22  illustrates  this 
occuirence.  The -9  to-  12  m/s  velcKities  in  the  ASR  vekx-'ity  field  (upper  right)  within  the  20.7  m/s 
alarm  were  used  by  the  detection  algorithm  even  though  they  conesponded  to  relatively  low  reflec¬ 
tivity  values  (upper  left).  The  nearest  con  esponding  TDWR  data  (bottom  panels)  and  UND  data  (not 
shown)  confirmed  microburst-strength  shear  was  not  present  in  the  alarm  region.  (Comparison  of 
the  ASR  and  TDWR  reflectivities  suggests  that  elevated-core  syndrome  may  also  have  contributed 
to  the  ASR  surface  velocity  overestimates  in  the  positive  velocity  portion  of  this  alarm.)  Possible 
fixes  for  this  problem  include  improving  the  ASR  Dual-beam  Processing  Algorithm,  or  rejecting 
velocity  estimates  associated  with  low-reflectivity  values  (e.g.  less  than  20  dBZ). 

Strength  and  range-dependent  area  thresholding  parameters  employed  by  the  microburst  algo¬ 
rithm  resulted  in  occasional  microburst  mis.ses,  mainly  early  in  the  lifetimes  of  microbursts  when 
they  were  quite  small  u[X)n  initial  impact  with  the  ground.  An  example  of  this  occurrence  is  provided 
in  Figure  23.  A  true  microburst  signature  was  evident  in  the  ASR  data  (top  panels)  within  the  17  m/s 
truth  polygon,  but  the  detection  algorithm  did  not  issue  an  alarm  because  the  area  of  the  “cluster”  of 
divergence  segments  was  slightly  less  than  the  required  1 .5  square  kilometers  for  an  event  at  that 
range  and  of  this  intensity.  An  alann  was  issued  by  the  detection  algorithm  a  few  seconds  later.  The 
conesponding  TDWR  tnath  data  is  shown  in  the  bottom  panels.  Short  delays  in  microburst  detection 
were  also  caused  by  cluster  age  requirements  imposed  by  the  microburst  detection  algorithm,  as 
alarms  were  not  declared  until  a  cluster  was  valid  for  three  consecutive  scans.  Conversely,  false 
alarms  occasionally  resulted  from  AMDA’s  coasting  of  alarms  after  a  cluster  disappeared. 

Microburst  asymmetry  accounted  for  some  of  the  false  alarms  and  microburst  misses  tabulated 
in  the  ASR-WSP  performance  scoring.  Occasionally  a  microburst  was  visible  in  the  ASR  data,  but 
not  in  the  TDWR  data,  when  the  ASR  had  a  more  favorable  viewing  angle  for  an  asymmetric  event. 
These  events  were  inappropriately  scored  as  false  alarms  because  single-Doppler  (TDWR)  data  was 
used  as  truth.  Figure  24  depicts  this  occurrence.  Both  ASR  (upper  panels)  and  UND  data  (often  used 
to  ascertain  presence  of  microburst  asymmetry)  in  the  lower  left  panel  indicated  microburst- 
strength  divergence  existed  in  the  1 8.7  m/s  alarm  region,  whereas  the  TDWR  (lower  right)  observed 
only  9  m/s  of  shear.  Viewing  angles  of  this  event  were  more  similar  between  the  ASR  and  UND,  so 
microburst  asymmetry  was  assumed.  Similarly,  when  a  microburst  was  not  evident  in  the  ASR  data, 
but  was  present  in  the  TDWR  data  because  of  microburst  asymmetry,  the  ASR-WSP  system  was 
penalized  with  misses.  The  lack  of  a  microburst-strength  signature  in  the  UND  velocity  field  con¬ 
firmed  the  asymmetry. 
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core.  Region  of  6  mis  approaching  veloci- 


Figure  21 .  Synthesized  RHI  scan  through  location  of  false  alarm  in  Figure  20.  Negative  velocities  associated 
with  elevated  reflectivity  core  dominate  ASR  surface  velocity  estimate  because  they  coincide  with  strongest 
signal  return. 
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-9  Microburst  Detection  Algorithm’s  use  of  velocity  estimates  in  areas  of  low  signal  return  caused  20  m/s  false  alarm.  TDWR  data  (bottom 
)  data  (not  shown)  confirmed  the  lack  of  microburst-strength  divergence  in  the  alarm  region. 


tiiiu*:  ‘‘ir7n7/n7  1  t  nm.'; 


Figure  23.  Missed  microburst  due  to  microburst  algorithm  area  thresholding.  Area  of  cluster  of  divergence  segments  in  ASR  data  (top  panels)  did  not  exceed 
minimum  required  for  event  at  this  range  and  of  this  intensity.  Nearest  TDWR  data  appears  in  bottom  panels. 
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The  main  types  and  eauses  of  AMDA  failure  have  been  described  above.  As  indicated,  many  ot 
the  observed  shoncomings  were  induced  by  physical  characteristics  ot  the  ASR  and  cunent  signal 
processing  techniques.  Improved  algorithm  pertbniiance  is  expected  after  enhancements  to  the  lat¬ 
ter.  Slight  modifications  to  the  algorithm  itself  might  also  augment  detection  performance. 


-49- 


6.  ONGOING  ALGORITHM  DEVELOPMENT 


There  are  two  sepitrate  paths  for  ongoing  ASR-9  microburst  algorithm  development.  The  first 
involves  relatively  minor  changes  to  the  existing  version  of  AMDA  to  improve  its  performance  us¬ 
ing  techniques  that  have  become  apparent  during  the  past  year  of  data  analysis.  The  second  involves 
alternative  algorithm  approaches,  motivated  by  recent  work  on  both  the  ASR-9  Gust  Front  Algo¬ 
rithm  [Delanoy,  Troxel,  1993]  and  the  Integrated  Terminal  Weather  System  (ITWS)  Microburst 
Detection  Algorithm  (Campbell,  1992], 

6.1  AMDA  ENHANCEMENTS 

AMDA  has  been  relatively  stable  since  August  1991,  but  no  relatively  complex  algorithm  is  ever 
completely  “finished”;  AMDA  is  no  exception.  As  the  volume  of  analyzed  data  has  continued  to 
grow,  problem  areas  have  been  uncovered,  and  minor  modifications  to  the  algorithm  are  continually 
being  evaluated  as  possible  solutions.  As  indicated  in  the  performance  evaluation  section,  noi.sy  ve¬ 
locity  estimates  in  areas  of  low  reflectivity  (<20  dBZ)  can  pose  a  problem.  The  noisy  estimates  are 
not  usually  coherent  enough  to  cause  a  false  alarm  on  their  own,  but  they  can  cause  a  significant  en  or 
in  the  velocity  estimate  for  true  microbui*sts  or  wind  shears.  The  most  common  result  is  that  an  event 
that  is  actually  only  a  wind  shear  alert  ( <  30  knots )  is  declared  a  microburst  alert  category,  contribut¬ 
ing  to  the  FAR  for  the  stronger  alamis.  A  modification  to  the  cun  ent  reflectivity  thresholding  logic  is 
being  evaluated  as  a  solution  to  this  problem.  A  threshold  would  be  applied  to  the  entire  velocity  field 
prior  to  the  segment  detection  process,  thereby  eliminating  the  possibility  that  valid  shear  segments 
contain  velocity  values  from  areas  of  questionable  signal  quality. 

A  second  problem  related  to  reflectivity  thresholding  is  the  performance  of  AMDA  in  a  dry  envi¬ 
ronment.  To  date,  all  ASR-9  microburst  data  have  been  collected  in  wet  environments.  This  will 
change  in  1993  with  the  relocation  of  the  mdar  to  Albuquerque,  New  Mexico.  To  detect  dry  micro¬ 
bursts  that  are  present  during  the  same  period  as  higher  reflectivity  storm  cells,  the  current  reflectiv¬ 
ity  thresholding  logic  will  need  further  modification.  Specifically,  to  avoid  the  removal  of  all  dry 
microburst  shear  segments,  the  reflectivity  threshold  parameter  will  have  to  be  set  on  a  storm  cell- 
by-storm  cell  basis  instead  of  utilizing  a  single  threshold  over  the  entire  scan. 

Though  the  final  microburst  outputs  are  being  temporally  smoothed,  the  temporal  stability  of  the 
alarms  can  sometimes  still  be  a  problem.  When  multiple  microbursts  are  present  and  in  close  prox¬ 
imity  to  one  another,  the  shapes  produced  by  the  current  algorithm  can  occasionally  change  signifi¬ 
cantly  from  scan  to  scan.  This  is  not  appropriate  output  for  a  rapid  update  data  stream  to  an  ATC 
environment,  and  becomes  even  more  of  an  issue  If  microburst  warnings  are  updated  more  rapidly 
than  the  30-second  interval  used  currently  during  operational  demos.  A  scan-to-scan  microburst 
shape  interpolation  strategy,  as  well  as  improvements  to  cunent  temporal  smoothing  logic  and  shape 
generation  algorithm  [Wilson,  et.  al.  1992],  are  being  studied  as  solutions  to  this  problem. 
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6.2  ALTERNATIVE  ALGORITHM  APPROACHES 


When  the  original  design  goals  of  AMDA  are  reexamined,  it  becomes  apparent  that  some  of 
them  are  no  longer  as  important  as  they  were  in  1988.  Changes  in  computer  hardware  and  software, 
as  well  as  in  the  definition  of  what  constitutes  the  "danger  zone’  of  a  microburst,  have  (xcuned  since 
that  time,  and  have  prompted  a  renewed  intere.st  in  alternative  approaches.  This  section  provides  a 
brief  histoi^  of  the  development  of  the  cunent  algorithm,  and  describes  some  perceived  shortcom¬ 
ings  of  AMDA  that  could  be  rectified  by  replacing  the  algorithm  with  caie  that  uses  more  refined 
image  pixxessing  techniques. 

6.2.1  AMDA  DEVELOPMENT  HISTORY 

When  the  TDWR  microburst  algorithm  was  first  developed,  it  was  decided  that  the  basic  goal 
was  todetect  the  peak  approaching  and  receding  velocities,  and  calculate  a  shear  value  (deltav)  as  the 
difference  between  the  two  peak  values.  A  second  design  consideration  was  the  decision  to  run  the 
algorithm  on  an  essentially  raw  (unsmoothed)  velocity  field,  again  in  the  interest  of  detecting  the 
absolute  peak  velocity  values.  The  combination  of  these  two  factors  made  the  implementation  of  a 
strictly  shear  field-based  algorithm  difficult  because  the  small  window  size  required  for  detection  of 
absolute  velocity  extrema  was  also  sensitive  to  noi.se  in  the  raw  vckxity  field.  For  this  reason,  a  seg- 
mcnt-ba.sed  point-to-point  search  algorithm  [Merritt,  1 987]  was  developed  instead  of  a  conceptual¬ 
ly  simpler  shear-field  based  algorithm  [Noyes,  1990],  and  while  the  resulting  algorithm  is  some¬ 
what  tricky  to  tunc,  it  has  been  very  successful  in  meeting  its  design  goal. 

The  ASR-9  microburst  detection  algorithm  was  not  required  to  mimic  the  TDWR  approach.  In 
fact,  the  design  considerations  were  slightly  different.  Due  to  spatial  smoothing,  the  ASR-9  velocity 
field  was  smoother  than  the  corresponding  TDWR  field,  and  an  initial  implementation  of  a  shear 
field-based  algorithm  was  successful  in  detecting  microbursts  of  vaiying  sizes  [Noyes,  1990J.  This 
algorithm,  however,  was  still  incomplete  in  early  1989,  and  a  requirement  existed  to  have  a  micro¬ 
burst  algorithm  running  in  the  WSP  that  summer.  To  accomplish  this,  a  simple  adaptation  of  the 
TDWR  algorithm  (AMDA)  was  developed  for  the  WSP.  During  early  1990  the  shear-based  algo¬ 
rithm  was  completed,  and  although  it  was  conceptually  simpler,  tests  conducted  using  recorded 
ASR-9  data  indicated  it  had  approximately  equivalent  performance  to  AMDA  and  was  more  com¬ 
putationally  intensive,  requiring  additional  processing  power  to  keep  pace  with  the  12  scan/min  an¬ 
tenna  rotation  rate.  Because  the  performance  of  the  two  algorithms  was  equivalent  and  AMDA  was 
already  running  in  real  time  in  the  WSP,  AMDA  remained  the  algorithm  of  choice  for  the  1990-1992 
time  period. 

6.2.2  MOTIVATION  FOR  ALTERNATIVE  ALGORITHMS 

Although  AMDA  runs  efficiently  and  has  performed  well  during  the  operational  test  periods, 
there  are  several  arguments  for  converting  the  algorithm  to  a  shear  field-based  approach.  The  main 
argument  is  that  the  existing  segment-based  ASR-9  and  TDWR  algorithms  do  not  necessarily  pro¬ 
vide  an  accurate  estimate  of  the  windspeed  loss  an  aircraft  will  experience  when  passing  through  a 
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microburst.  In  many  cases  durinii  the  operational  demonstrations,  pilots  flying  through  wind  shear 
events  would  not  encounter  the  predicted  loss  value  because  the  cunent  warning  supplies  only  the 
overall  pcak-to-peak  vekxity  diflerential,  with  no  information  about  the  rate  of  change.  For  exam¬ 
ple,  the  loss  value  provided  t(^  pilots  for  a  4()-knot  microburst  with  a  diameter  of  4  km  is  identical  to 
the  loss  for  a  4()-knot  niicroburst  with  only  a  1  km  diameter,  but  an  aircraft  will  be  affected  quite 
differently  in  the  two  cases.  Funhermorc,  two  microbursts  with  identical  size  and  peak-to-peak  ve¬ 
locity  difference  may  also  be  quite  different  if  the  shear  is  spread  out  fairly  evenly  in  one  case  but  not 
in  the  other.  A  recent  study  conducted  foi  TDWR  compared  the  issued  microburst  wiunings  with  the 
actual  performance  degradiUion  encountered  by  a  research  aircraft  flying  through  over  60  micro¬ 
bursts  in  Orlando.  Florida  | Campbell.  1 992 1.  The  performance  degradation  on  board  the  aircraft  was 
measured  using  a  hazard  index  known  as  "F-factor”  (Bowles,  1990],  which  takes  into  account  the 
overall  energy  balance  of  the  individual  aircraft  as  well  as  the  sunounding  environment.  To  make  the 
ground-based  radar  warnings  compare  favorably  with  the  measured  F-factor  data,  it  was  ncccssai7 
to  compute  an  equivalent  F-factor  using  a  shear  field.  In  light  of  this  result,  it  would  be  desirable  to 
have  the  sheai"  field  available  to  the  algorithm,  so  that  any  future  upgrades  to  the  algorithm  outputs 
could  include  the  shetu-  information. 

Recent  work  in  the  area  of  gust  front  detection  has  resulted  in  the  new  ASR-9  Maehine  Intelli¬ 
gent  Gust  Front  Algorithm  (MIGFA)  [Delanoy,  Troxel,  1993],  This  algorithm  has  its  roots  in  ad- 
vaneed  image  processing  techniques  (also  utilizing  a  shear  field)  and  is  a  dramatic  improvement, 
both  in  performance  and  in  ea.se  of  modification,  over  the  previous  segment-based  algorithm.  The 
major  challenge  when  detecting  gust  fronts  with  an  ASR-9  is  the  ability  to  extract  gust  front  thinlines 
from  low  SNR  radar  data.  This  is  similar  to  the  problem  of  "dry”  microburst  detection  with  an 
ASR-9,  and  given  the  success  of  MIGFA  in  this  regard,  it  may  prove  necessaiy  to  test  this  approach 
for  dry  microburst  detection  as  well.  Since  the  WSP  will  already  contain  one  algorithm  of  this  type, 
this  strategy  would  also  make  the  overall  architecture  of  the  WSP  more  consistent,  and  therefore 
easier  to  produce  and  maintain. 

No  major  eommitment  has  yet  been  made  to  replaee  AMDA  with  a  more  flexible  and  eomputa- 
tionally  intensive  shear  field-based  algorithm.  Much  of  the  background  work  has  already  been  per¬ 
formed  however,  and  the  effort  required  to  finish  the  task  is  certainly  not  overwhelming.  The  current 
generation  of  processing  hardware  provides  sufficient  performance  to  do  the  job,  and  a  new  “ma¬ 
chine-intelligent”  shear-based  algorithm  could  be  ported  to  the  WSP  given  sufficient  resources. 
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7.  SUMMARY 


AM  DA  evolved  from  the  TDWR  microburst  detection  algorithm  and  has  been  optimized  tor  the 
ASR-9  radar.  Like  the  TDWR  algorithm,  AMDA  detects  shear  .segments,  and  groups  the  segments 
into  wind  shear  alens  and  microburst  ahu'ms.  Unlike  the  TDWR  algorithm,  AMDA  utilizes  a  coordi¬ 
nate  system  shift  to  process  oveihead  data  in  a  convenient  manner,  and  also  makes  use  of  the  retlec- 
fivity  field  to  help  eliminate  fal.se  ahums  caused  by  precipitation  aloft.  The  algorithm  is  efficient, 
cuiTcntly  running  in  real  time  on  a  5  MIP,  68()3()-bascd  single  board  computer  with  a  worst-case 
CPU  usage  of  approximately  30  percent. 

AMDA  has  been  successful  during  the  operational  evaluation  test  peritxls  in  Orlando  during  the 
past  two  summers,  with  an  average  detection  petformance  of  97  percent  and  a  falsc-alarm  rate  not 
exceeding  1 0  percent  for  events  stronger  than  30  knots  and  within  8km.  AMDA  has  not  yet  under- 
g(Mie  testing  in  a  dry  environment,  and  changes  to  the  algorithm,  specifically  in  the  reflectivity 
thresholding  area,  may  be  necessary  to  achieve  acceptable  peifomiance  in  this  case.  The  move  to 
Albut|ucrquc  in  1993  should  provide  the  neces.saiy  diy  microburst  data  to  evaluate  the  algorithm  in 
this  regard. 

Increased  understanding  of  the  operational  issues  associated  with  microburst  detection,  and  the 
successful  implementation  of  the  Machine  Intelligent  Gust  Front  Algorithm  have  prompted  a  re¬ 
newed  interest  in  alternative  algorithm  approaches.  With  the  advent  of  faster  computer  hardware, 
real-time  implementation  of  these  computationally  intensive  approaches  is  now  becoming  more 
practical,  and  it  is  worth  evaluating  these  methods  when  performing  further  algorithm  development. 
Potential  benefits  of  a  “Machine  Intelligent”  microburst  algorithm  include  improved  performance, 
particularly  in  challenging  meteorological  environments,  and  greater  consistency  among  the  WSP 
mici  oburst  and  gust-front  algorithms. 
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Appendix  A  AMDA  Parameters 


This  appendix  contains  a  list  of  all  the  adjustable  parameters  used  by  AMDA.  The  current  settings  are  shown 
in  parentheses  to  the  right  of  each  parameter. 


Segment  Detection  Parameters 

SEG_MAX_1ST_INC  (5.0  m/s) 

SEG_MIN_LENGTH  (.96  km) 
SEG_MIN_SHEAR  (7.8  m/s) 

MAX_SUM_DECR£ASES  (5.0  m/s) 

SUM_DECREASES_ADJ  (1.0  m/s) 


SEG_MAX_BAD_VALS  (2) 
SEG_MIN_SLOPE  (1.8  (m/s)dan) 


Maximum  velocity  differential  between  first  two  gates  of 
a  segment. 

Minimum  length  for  a  segment  to  be  considered  valid. 

Minimum  velocity  differential  for  a  segment  to  be  consid¬ 
ered  valid. 

Maximum  value  for  sum_decreases.  Values  exceeding 
this  will  cause  the  segment  growing  process  to  terminate. 

Adjustment  value  to  apply  lo sum_decreases  when  a  veloc¬ 
ity  increase  is  encountered  and  sum_decreases  is  non¬ 
zero. 

Maximum  number  of  bad  values  allowed  in  a  segment. 

Minimum  slope  requirement  for  segments.  Segment  re¬ 
gions  with  slope  less  than  this  are  trimmed  away. 


Segment  Temporal  Smoothing  Parameters 

SEG_AGE_MIN  (3  scans)  Minimum  value  for  the  segjige  variable  before  a  segment 

can  be  inserted  into  a  scan  where  it  is  missing. 

SEG_AGE_MAX  (7  scans)  Maximum  value  allowed  for  the  segjige  variable.  The 

formula  (SEG_AGE_MAX  -  SEG_AGE_MIN)  repre¬ 
sents  the  maximum  number  of  segment  insertions  that  can 
be  done  after  a  segment  disappears. 

Segment  Association  Parameters 

SEG_OVERLAP_PERCENTAGE  (50%)  Percentage  of  the  segment  to  use  in  the  overlap  test  In  this 

case,  only  the  middle  50%  of  the  segment  is  used  for  over¬ 
ly  testing. 

CLU_MIN_SEGS  (4)  Minimum  number  of  segments  in  a  cluster. 

CX.U_MIN_AZ_WINDOW  (3)  Minimum  azimuth  window  to  use  during  segment  associa- 

tiwi.  The  azimuth  window  varies  with  range,  and  this 
group  of  4  parameters  defines  the  endpoints  of  a  linear  azi¬ 
muth  window  vs.  range  equation. 

CLU_MIN_AZ_WINDOW_RANGE(  10.0km)  Range  at  which  to  use  the  minimum  azimuth  v/indow. 

CLU_MAX_AZ_WINDOW  (11)  Maximum  azimuth  window  to  use  during  segment 

association. 
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CLU_MAX_AZ_WlNDOW_RANGE(O.Okm)  Range  at  which  to  use  the  maximum  azimuth  window. 


Reflectivity  Thresholding  Pat-ometers 

REFL_THRESH_CORE_DBZ  (40.0  dBZ) 
REFL_THRESH_SEG_LEN_1N_CORE(50%) 

REFL_THRESH_DIST  (0.5  km) 


Segment  Density  Thresholding  Parameters 

SEG_DENSITY_THRESH  (0.5)  Segment  density  thresholding  parameter.  Ouster  seg¬ 

ments  in  regions  where  the  segment  density  is  below  this 
fraction  are  discarded. 

Area  Thresholding  Parameters 

CLU_AREA_MIN  (0.5  km^) 

CLU_AREA_RANGE_MIN  (2.0  km^) 

CLU_AREA_RANGE_MAX  (16.0  km2) 

CLU_AREA_RANGE_SLOPE  (0.1  km^/km) 

CLU_AREA_SHEAR_A  (30) 

CLU_AREA_SHEAR_B  (2.5) 

CXU_AREA_SHEAR_C  (-1.0) 

Cluster  Temporal  Smoothing  Parameters 

CLU_MIN_AGE  (3  scans)  Minimum  cluster  age.  Younger  clusters  are  not  ouQtut  as 

microburst  events. 

CLU_COAST_SCANS  ( 12  scans)  Number  of  scans  to  coast  a  cluster  if  it  is  dropped. 


Baseline  value  of  area  threshold. 

Minimum  range  fw  area  threshold  range  correction. 
Maximum  range  for  area  threshold  range  correction. 
Rate  of  change  of  area  threshold  with  respect  to  range. 
A  Coefficient  for  area  threshold  shear  ejection. 

B  Coefficient  for  area  threshold  shear  correction. 

C  Coefficient  for  area  threshold  shear  correction. 


Core  reflectivity  used  for  reflectivity  threshold. 

Percentage  of  segment  length  that  must  lie  within  a  region 
of  core  reflectivity  in  order  for  the  segment  to  be  consid¬ 
ered  valid. 

Distance  used  in  the  last  stage  of  the  reflectivity  threshold¬ 
ing  process.  Segments  within  this  distance  from  a  segment 
that  has  passed  the  reflectivity  test  are  also  considered  val¬ 
id. 
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